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Electronic properties of a graphene antidot in magnetic fields 
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We report on several unusual properties of a graphene antidot created by a piecewise constant 
potential in a magnetic field. We find that the total probability of finding the electron in the 
barrier can be nearly one while it is almost zero outside the barrier. In addition, for each electron 
state of a graphene antidot there is a dot state with exactly the same wavefunction but with a 
different energy. This symmetry is a consequence of Klein tunneling of Dirac electrons. Moreover, 
in zigzag nanoribbons we find strong coupling between some antidot states and zigzag edge states. 
Experimental tests of these effects are proposed. 

PACS numbers: 



I. INTRODUCTION 

Bul k La ndau energy levels of graphene, i?„ = 
^Ec^y2\n\, have different magnetic field B and quan- 
tum number n dependence than those of ordinary two- 
dimensional Landau levels, = hLdf.{n+ 1/2) (Here the 
cyclotron energy hujc oc B, where B is the magnetic field. 
The characteristic energy of Dirac electrons Ec cx VB). 
The wavefunctions of these two systems have, respec- 
tively, two-component and one-component structures 
In the presence of a dot or antidot potential these states 
of the two systems will be perturbed rather differently. 

Electron wavefunctions of a graphene quantum dot 
have been investi gate d actively both with and without 
a magnetic field In absence of a magnetic field an 
electron cannot be localized due to Klein tunneling and 
only quasibound states are allowed. However, a perpen- 
dicular magnetic field enhances the localization of the 
wavefunctions and both bound states and quasibound 
states occur. 

Recently properties of antidot lattices [T3l - [T9| have 
been explored actively. But properties of a single an- 
tidot in a magnetic field has not been studied thoroughly 
yet. We believe that an antidot created by cutting a hole 
out of the graphene sheet [H, [1^1 and an antidot induced 
by an electrostatic gate can have different properties. In 
an antidot created by holes skipping orbits that encircle 
an integer number of elementary quantum flux defined by 
the antidot are allowed [2C|: BA ^ mipQ, where A is the 
enclosing area of the antidot, m is the orbital quantum 
number and 0o is the magnetic fiux quantum. However, 
some graphene states of a gate induced antidot potential 
may not obey this quantization rule since Klein tunnel- 
ing allows significant penetration into the barrier. Hence, 
graphene antidots may not always support skipping or- 
bits at the edge of an antidot. Also we expect an antidot 
defined by a piecewise constant potential may have sev- 
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FIG. 1: Schematic drawings of some typical probability wave- 
functions of antidots and dots. We consider piecewise con- 
stant antidot and dot potentials Vair) and Vdij) satisfying 
Va{r) + Vd{r) = V. (a) and (b) display dot and antidot 
states with the same wavefunction but with different ener- 
gies Ed ~ V — Ea and Ea- Solid (dashed) line represents the 
probability wavefunction of A (B) component. Naively one 
may expect that the energy of the antidot state of (a) should 
be larger than the dot state by amount V since the antidot 
state is in the barrier with the height V , but this reasoning is 
actually incorrect, (c) and (d) are same as in (a) and (b) but 
with Va{r) + Vd{r)=0. Antidot and dot states have energies 
are Ea and Ed = —Ea- (e) and (f) are same as (a) and (b) 
but wavefunctions are localized outside the antidot and dot 
regions, (g) shows tight-binding probability wavefunction of 
an antidot state of a zigzag nanoribbon. 



eral interesting properties in the presence of a magnetic 
field. This is because of the following feature of the Dirac 
equation: eigenstate wavefunctions of the Dirac equation 
are also eigenstates of Schrodinger-like equation, where 
constant potentials V and eigenvalue E appear together 
in the effective energy as E'^^ = {E — V)'^. This is an in- 
teresting feature of Dirac electrons but its consequences 
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have not been fully explored in a magnetic field. Nature 
of antidot states of a nanoribbon with zigzag edges may 
also be interesting. A nanoribbon has surface edges with 
various localization lengths j2ll[2^. It is unclear whether 
these chiral zigzag edge states may be coupled to confined 
non-chiral antidot states. 

In this paper we report on our investigation of these 
issues. We consider cylindrically symmetric and piece- 
wise constant antidot potentials. We find that for each 
electron state of a graphene antidot there is a dot state 
with exactly the same wavefunction but with a differ- 
ent energy, see FigslIJa)-(f). This symmetry is a con- 
sequence of the appearance of E^g in the Schrodinger- 
like equation. We find that the eigenstates may be di- 
vided into three classes: those that are localized inside 
the antidot(FiglTIa)), those that have significant weights 
at the boundary (Fig [TJc)), and those that are localized 
outside(Fig[He)). Interestingly we find that probability 
of finding the electron in the classically forbidden region 
of the barrier can be almost one, as shown in Figs[TJa) 
and (f ) . It is a counter example to the usual expectation 
that the probability wavefunction is larger in a potential 
well than in a barrier. This is a consequence of an in- 
terplay between strong Klein tunneling and localization 
of Landau level wavefunctions. In zigzag nanoribbons 
we find strong coupling between some antidot states and 
zigzag edge states, see FiglTJg). This effect is unique to 
antidots of zigzag nanoribbons. Experimental tests of 
these effects are proposed. 



II. SYMMETRY BETWEEN ANTIDOT AND 
DOT STATES: 

We consider a cylindrically symmetric and piecewise 
constant potential of an antidot or dot: where the radius 
of the antidot and dot is R: V{r) = Vi for r < R and 
V{r) = Vii for r > R. Around the K point of the Bril- 
louin zone the Dirac equation of an antidot or dot has 
the form 



H = vpa ■ {p+ -A) + V{r) 
c 



(1) 



with the elementary charge e > 0, the Fermi velocity wf, 
the Pauli spin matrices a = (ca;, ctj,, CTz), and magnetic 
vector potential A = ^{—y,x,Q). Magnetic field B is 
along z-axis and is perpendicular to the graphene layer. 
Since the Hamiltonian commutes with the angular mo- 
mentum operator = —id^ -\- <Tz/2 the eigenstates can 
be written in polar coordinates as 



XA{r) 



(2) 



For each angular momentum number j = 
±1/2, ±3/2, ±5/2, ••• there are numerous excited 
eigenvalues Ej„i, labeled, in increasing order, by 
TO = • • • , -2, -1, 0, 1, 2, • • • {Ej^„, > for TO > and 



It follows from Eq. ([T]) that A and B components of the 
wavefunction "00-, a with energy Ea satisfy 



i'a.c = {{Ea - 14)' T £^e)V'.- (3) 



I, II stands for r < R or r > R 
denotes for A (B) component. 

and £ — ^ / ttht stand for the 



where the subscript a = 
and the index a = —{+] 

The constants E^ ^ ^ 

characteristic energy scale of Dirac electrons in magnetic 
fields and the magnetic length. If the potential is not 
piecewise constant another term containing ^ will be 
present in this equation. The equation has the same 
mathematical structure as the Schrodinger equation of 
the two-dimensional Landau levels, but there is an im- 
portant difference: the eigenenergy Ea and the constant 
potential 14 appears together as E^^^ = [Ea — 14)'- For 
V{r) = we obtain the usual result of two-dimensional 
Dirac electrons in a magnetic field, E^ — 2E^n. Now 
consider a cylindrical antidot and dot with piecewise con- 
stant potentials Va{r) and Vd{r) such that the sum of 
them is a constant Vd{r) + Va{r) — V. Then, it can be 
shown from Eq.Q that, when the antidot potential Va{r) 
has an eigenenergy Ea, the dot potential Vd{r) will have 
an eigenenergy Ed — V — Ea with the same eigenstate. 
Several examples for different pairs of dot and antidot 
are shown in FiglTJ 
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FIG. 2: Energy spectrum of antidot as a function of magnetic 
field. 




FIG. 3: Plot of the probability wavefunctions of state j — 5/2 
at two different value of B. Inset: its energy level decreases 
as B increases. 



Ej^rn < for TO < 0). 
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FIG. 4: Probability wavefunctions of antidot with radius R at 
B = 10.28T(£ = 80A), A and B components are represented 
by solid and dashed lines. 
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FIG. 5: Same as in FigHbut for j = ±3/2. 



III. ANTIDOT EIGENVALUES AND 
EIGENSTATES 

The properties of probability wavefunctions and 
eigenenergies of a graphene antidot are rather different 
from those of the ordinary antidot. This is because skip- 
ping orbits are not well defined in graphene antidot due 
to strong Klein tunneling in the presence of a magnetic 
field. The eigenstatcs and eigenvalues of the antidot may 
be found by matching solutions inside and outside of the 
antidot at the boundary r — R: the exact eigenstate 
wavefunctions of Eq. ([T]) for a cylindrical potential can be 
described by confluent hypergeometric functions 0, 0| h> 
stead of Laguerre polynomials of bulk Landau levels. The 
energy spectrum is displayed in Figl2] for Vi = 0.26eV, 
Vii = 0, and R = 200A. In the strong magnetic field 
limit R/£ 00 we have tested numerically that the bulk 
graphene Landau level energies are recovered. 

Many positive (negative) energy levels in Figl2]increase 



(decrease) as vi3. Their wavefunctions are localized out- 
side the antidot and are nearly unaffected by the antidot 
potential. An example is shown in Figdje) . Other energy 
levels display different dependence on B. Their wave- 
functions are affected by the antidot potential and have 
significant weights on the boundary r — R. An example 
is given in FigEl As one can expect from the usual Lan- 
dau level physics, the effective energy Ecs = E—V of this 
level increases as B increases, and its wavef unction be- 
comes more confined in the antidot region. However, this 
implies that the cigenenergy E decreases with increasing 
B, as shown in the inset of Fig [31 

We have investigated the eigenfunctions for different 
values of E and j. The first three lowest positive energy 
states for several values of j are shown in Figs|4] and [S] 
In contrast to ordinary tunneling physics, we observe a 
significant penetration into the antidot region: note that 
for j = —1/2 and —3/2 the probability wavefunctions of 
A component are smaller than probability wavefunctions 
of B component in the antidot region. For j = 1/2 and 
3/2 the opposite is true. The probability wavefunctions 
of the lowest energy states with j = 1/2 are plotted in 
FigHl (d): in the limit r — )■ the wavefunction XAir) 

is proportional to while XBir) is proportional to 

re~^ . However, for j = —1/2 the opposite is true, as 

^2 

shown in Fig|4](a): XAir) is proportional to re~«^ and 

^2 

Xsir) is proportional to e~«^. We find that as \i\ in- 
creases the mean radii of the wavefunctions ^ jm increase. 



When the mean radius 



is greater than R the 



probability wavefunctions will be peaked in the barrier. 
In a quantum well the wavefunction for a large for 
example, —51/2, is thus strongly localized in the harrier, 
as shown schematically in Figlljf). This unusual effect 
is possible because the constant potential of the barrier 
appears only in [E — I/)^(see Eq.Q). The wavefunction 
of this state has a rather different form in comparison 
with quasibound states at zero rnagnetic field with a sig- 
nificant weight inside the dot [2!. [lOl [T2| . 




FIG. 6: Occupation probabilities on A and B carbon atoms 
of a state near the Dirac point with E = — O.OleV. Small 
black dots indicate lattice sites and red (blue) dots indicate 
values of probability wavefunction on A (B) carbon atoms. 
The circle represents the antidot. The unit of x and y axis is 
a = 2.461. 
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(a) 



(c) 



(b) 





FIG. 7: (a) and (b) display probability wavefunction of a state 
with i5 = 0.18 eV and its pseudospin picture, (c) and (d) show 
probability wavefunction of a state with E — —0.04 eV and 
its pseudospin picture. 



IV. ANTIDOT STATES OF ZIGZAG 
NANORIBBONS 



the corresponding result of Dirac electrons, see FigslSfa) 
and (d). Tight-binding and Dirac electron models thus 
yield qualitatively similar results. 



V. DISCUSSIONS AND CONCLUSIONS 




In zigzag nanoribbons it would be interesting to in- 
vestigate the interplay between Klein tunneling and 
backscattering [25l - l27| . We adopt the following model for 
a nanoribbon with an antidot. The zigzag edges are hor- 
izontal along the x-axis and periodic boundary condition 
is imposed on the left and right edges. The horizontal 
and vertical lengths of the nanoribbon are — 130. 35A 
and Ly = 130.641. A magnetic field of B = 20T is ap- 
plied perpendicular to the graphene sheet. The radius 
of the antidot is 19.68A and its the potential height is 
V = leV. We solve for eigenstates and eigenvalues in 
the usual tight-binding Hamiltonian. FiglH] displays the 
occupation probabilities on A and B carbon atoms of a 
eigenstate close to the Dirac point with nearly zero en- 
ergy (its probability wavefunction is shown in Figlljg)). 
This state is a mixture of zigzag edge and antidot states. 
Opposite chiralities are found on the opposite sides of the 
zigzag edges. Since occupation is chiral along the zigzag 
edges no current flows along the edges. This mixed state 
is thus not a magnetic edge state. In the absence of an 
antidot some states near the Dirac point can have long lo- 
calization l eng ths while states away from the Dirac point 
do not [111, [l^- Some of these states with localization 
lengths comparable to the width of the nanoribbon cou- 
ple with the antidot, and significant probability wave- 
function can be present both at the edges and near the 
antidot. 

On the other hand, probability wavefunctions of eigen- 
states away from the Dirac point are not peaked along the 
edges. Their properties are displayed in Figl?! FigsITJa) 
and (b) show an antidot state with the probability wave- 
functions that are peaked at the origin of the anti- 
dot. Note that a similar feature is observed with Dirac 
electrons, as shown in Figs|31Ja) and (d). In contrast, 
FigslTJc) and (d) display a state with the probability 
wavefunction that is zero at the origin, consistent with 



FIG. 8: (a) and (b) display probability wavefunctions of a 
degenerate pair of K and K' with E — 0.188eV. A component 
is dashed line and B component is solid line, These probability 
wavefunctions are obtained solving Dirac equations. 



0.004 




0.000 1 ■ ' ■ ' ■ ' 1 

12 3 

r/R 

FIG. 9: The solid line shows the total probability wavefunc- 
tion at lattice sites with E = 0.188eV obtained by solving 
Dirac equation and the dashed line displays the total prob- 
ability wavefunction with E = 0.175eV obtained by solving 
the tight-binding Hamiltonian. 



In out treatment of Dirac equations we ignored the cou- 
pling of states of K and K' points due to the presence of 
a sharp boundary[2^ in the antidot potential. In the ab- 
sence of the sharp boundary states of K and K' points are 
degenerate with wavefunctions of A and B components 
exchanged, as shown in FiglS] (The radius of the antidot 
is 19.68A and its potential height is y = leV with the 
magnetic field B = 20T). The coupling between K and 
K' points can be included in a tight-binding model, and 
our calculation shows that this degeneracy of Dirac equa- 
tions is broken slightly: for example, E = 0.188eV splits 
into E — 0.175eV and O.lSleV, Moreover, the wave- 
functions of tight-binding and Dirac equation are rather 
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close to each other, see Fig|9l Furthermore, inside the 
antidot A and B components of nearly degenerate tight- 
binding wavefunctions are approximately similar to those 
of degenerate solutions of Dirac equations, see Figsl5]and 
[TOl A components are dominant in FigsHl^a) and [TOla) 
while B components are dominant in FigsEl^b) and iroT b). 
These tight-binding results demonstrate that the neglect 
of the coupling between K and K' points in the approach 
of Dirac equations is a reasonable approximation. 



(a) 




In this paper we have discussed several unusual prop- 
erties of tunneling physics of step-like graphene poten- 
tial barriers in magnetic fields. Scanning tunneling 
microscope may be used to test strong Klein tun- 
neling effects in antidots since the electron density in 
the antidot should be non-zero in contrast to an ordi- 
nary antidot. Infrared optical transitions f2^ may be 
used to test the symmetry between antidot and dot wave- 
functions discussed in this paper: Consider an occupied 
level Ea,i and an unoccupied level i?a,2 in an antidot. 
Then, according to the proposed symmetry, the absorp- 
tion transition Ea^i — >■ £^a,2 of the antidot will have the 
same transition probability as the absorption transition 
Ed,2 = V" - £^0,2 Eds = V - Ea,i of the corresponding 
dot system. We have also found a significant penetration 
of wavefunctions deep into the barrier region, as shown 
in Fig. 1(f). This effect is absent in the usual quantum 
Hall physics in Hall bar geometry, and it may be worth- 
while to explore how it may affect quantum Hall edges 
of graphene. 
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FIG. 10: Occupation probabilities on A and B carbon 
atoms of a state near the Dirac point with E=0.181eV(a), 
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